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SIMULATION OF THE PRIMARY LOOP OF A NUCLEAR POWER PLANT WITH A SMALL, 

GENERAL PURPOSE ANALOG COMPUTER 

INTRODUCTION 

This study describes the simulation of the 4 'primary loop" of a nuclear 
power plant on a twenty amplifier, general purpose analog computer. This 
simulation, which normally requires many more amplifiers, was performed 
using a passive network to simulate reactor kinetics and a capacitor 
memory to simulate transport delay. These special components did not 
impair the versatility of the general purpose computer since they are inter- 
changeable with the computer's normal computing components. 

In this simulation, the primary loop was considered to be operating under 
steady- state conditions at one-half of its maximum power. The objects 
of the simulation were: to determine system response to demand power 
changes and step inputs in reactivity; to observe reactor response to con- 
trol system failure during a demand power change; to determine whether 
the reactor can be controlled manually; and to illustrate reactor behavior 
when 14 scram rods" are inserted in the system. The results are sufficiently 
accurate for an approximate analysis of the over- all behavior of the sys- 
tem and are useful in practice. 

The simulation and special purpose equipment will be described after 
considering first the physical system and the mathematical model in 
detail. In this way, a more complete explanation of their function can be 
given. This will be followed by a discussion of the mechanization of the 
model and the results of the simulation. 
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SIMULATION OF THE PRIMARY LOOP OF A NUCLEAR POWER PLANT WITH A SMALL, 

GENERAL PURPOSE ANALOG COMPUTER 



PHYSICAL SYSTEM 

The physical system to be simulated is shown 
schematically in Figure 1. The fission process, 
which takes place in the U 235 fuel elements, gen- 
erates a large amount of heat energy. This energy, 
which raises the temperature of the fuel elements, 
the moderator, etc., is transferred to a coolant 
flowing through the core of the reactor. This 
coolant, a pressurized liquid, leaves the reactor 
through an outlet plenum (mixing chamber) at an 
elevated temperature and is transported to a steam 
generator. In the steam generator, the coolant 
transfers heat to feedwater or low quality steam, 
converting it to high pressure steam, which flows 
in a secondary loop. This steam is then converted 
to electricity in the conventional manner or put 
to some other use. The coolant leaving the steam 
generator is transported back to the core of the re- 
actor, entering through an inlet plenum chamber. 



The control system of the reactor consists of a 
comparator, a controller, and a control rod drive 
unit. 

The comparator, which compares actual reactor 
power level with desired power level, generates 
an error signal corresponding with any detected 
difference. This error signal is fed to a controller 
which in turn positions the reactor-control-rods to 
correct the difference. The position (depth of in- 
sertion) of the neutron- absorbing control rods, 
which are located adjacent to the fuel elements, 
governs the number of fissionable neutrons pro- 
duced per fission. The neutron density in the fuel, 
in turn, governs the amount of energy released 
by the fission process. This energy release is 
the determining factor of the actual power level 
of the reactor. 
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Figure 1. Primary Loop of a Nuclear Power Plant, Simplified Schematic 



The purpose of the ' 6 scram' 9 system is to position 
the neutron- absorbing rods in the reactor rapidly 
when the multiplication of neutrons becomes un- 
controllable. These scram rods, which are spring- 
loaded or magnetically held, absorb enough neu- 
trons to stop the fission process safely. 



The reactivity of the system, 6 k(t) , is normally 
defined in terms of the multiplication factor, k(t) , 
by the equation 



Sk (t) = 



MtH 
k(t) 



(3) 



The various components of the 6 'primary loop", and 
their interdependence, will be discussed in the 
formulation of the mathematical model. 

Since this simulation dealt with the behavior of 
the loop for various disturbances, and not with 
equipment design, the mathematical model was de- 
rived in terms of average quantities. This per- 
mitted mathematical definition without the use of 
partial differential equations. 

In addition, certain restrictions were placed on 
the system in order to stay within the limits of a 
twenty-amplifier computer. First, coolant flow 
rate was assumed constant during all transient 
periods of operation. Second, average steam tem- 
perature in the generator was assumed constant 
at all times. This restriction limited the simula- 
tion to the primary loop. 



MATHEMATICAL MODEL 



wfrere k(t)is the number of fissionable neutrons 
produced per fission. In order for the chain re- 
action to be self-sustaining, this multiplication fac- 
tor must be unity and reactivity must be zero. When 
these conditions are met the reactor is said to be 
in a "critical" state, and will continue the chain 
reaction or fission process for long periods of 
time or until stopped. When reactivity is greater 
than zero [k(t) > l] , the reactor is in a " super- 
critical" (disastrous) state; when less than zero 
[k(t)<l] , the reactor is in a "subcritical" state 
and will eventually cease to produce neutrons and 
stop the process. 

For computational purposes, the system reactivity 
is normally expressed as the linear combination 
of its contributing reactivities (3, 7), namely: 



8k(t) = Sk f (t) + Sk (t) 



(4) 



+ 8k f (t) + 8k c (t ) 



The following areas of the primary loop must be 
defined mathematically to achieve a successful 
simulation: 

1„ Reactor Kinetics — The kinetic equations for 

nor - 1 

\J* 0 ° fuel can be found in several references (1, 2, 
3 , 4 , 6, 7 ). These equations, which take into ac- 
count delayed as well as prompt neutrons, are 
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(1) 



Since changes in 5k^(t) with time depend upon fuel 
depletion, breeding, etc., it can be considered as 
constant for short transient periods. This is also 
true for neutron- absorbing reactor poisons such 
as Xe 135 . Their concentration changes have a 
negligible effect on over-all reactivity compared 
to the effects of fuel temperature and control rod 
position. It should be noted, however, that the 
effect of poisons is quite important during simula- 
tions of start-up and shut-down procedures. 



The reactivity contribution made by fuel tempera- 
ture is represented mathematically by the Equa- 
tions (3, 6, 7) 
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In this study, the temperature coefficient of re- 
activity is negative (a fuel characteristic) 
and the system is self regulating, i.e., as the re- 
actor heats up, its reactivity is reduced. The con- 
tribution of reactor poisons and control rod posi- 
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tion are also negative, and tend to reduce system 
reactivity and produce a "subcritical" state. The 
built-in reactivity of the fuel has a positive effect 
on reactivity. 

Accordingly, Equation 4 may be rewritten in the 
form 



This transfer of energy is usually described by a 
partial differential equation; in the interest of 
simplicity, it was linearized to obtain Equation 10 
(the mathematics involved are shown in Appendix A) . 
A byproduct of this linearization is that inlet, 
average, and outlet temperatures of coolant flowing 
in core are related by the equation 



Sk(t) = k + QCT f (t) +,i(t) 



Here 

Mt) = 8k ( t)- 8k (0) 
c c 



(6) 



(7) 



and 



k = 8k c (0) + 8k p + 8k f - * T Q (8) 



Equation 6 is the most convenient form of the re- 
activity equation to mechanize. 

2o Fuel Element Heat Trans — Considering re- 
actor fuel elements, cladding, etc., as a whole, 
the differential equation governing heat transfer 
in the fuel elements is 
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In deriving Equation 9, the rate of heat generation 
in the fuel was assumed proportional to neutron 
density, n (t). The proportionality constant, AH^, is 
defined as the heat of fission. 

3 p Fuel Element, Coolant Heat Transfer — The 
equation describing transfer of heat energy be- 
tween fuel elements and flowing coolant can be 
given, verbally, as: total heat energy accumulated 
by coolant in reactor equals the difference between 
heat input by fuel and total heat removed by coolant, 
or 



dT.(t) 



M c C c^- = UA[T f (t)-T c (t)] 
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(10) 
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(11) 



This equation may be combined with Equation 10 to 
obtain 



M C 
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c c dt 
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(12) 



or 



dT (t) /UA+2W C \ 

+ ( .. - ° ) T c (t) = UAT f (t) 



dt 



M c c c 



(13) 



+ 2W c C c T ic (t) 



which is the equation for a first order lag. 

4„ Mixing in Plenum Chambers — These equa- 
tions were derived by treating the chambers as 
stirred tanks. Thus 



M 



and 



M 



dT|c(,, = W c [T i(t ,-T lc ,.,] W 



dt 



For convenience, Equation 14 can be combined 
with Equation 11. 



(16) 



M 0 — = W c [2T c (t)-T ic (t)-T 0 (t)J 



5o Transport Delay — Generation of this function 
is necessary to simulate flow in the piping sec- 
tions between outlet plenum chamber and inlet to 
steam generator, and output of steam generator 
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and inlet to plenum chamber. Coolant will be de- 
layed a given length of time, t Q9 as determined 
by its mean velocity, v , and length of the piping 
section, £ Q . The function which must be generated 
to satisfy this requirement is 



is determined by a control program. A typical 
control program is to maintain constant average 
temperature (4). This program can be achieved 
by means of a proportional plus integral controller 
described mathematically by 



and 



Here 



V*> = T o (t - V 



T i (t) = Vt-T,) 



(17) 



(18) 



T o= 4r and t. = 4 1 (19) 



6oHeat Transfer in Steam Generator — This 
equation is very similar to the one governing heat 
transfer between fuel and coolant in reactor core. 
That is, total heat energy accumulated by coolant 
in the steam generator equals difference between 
heat input and total heat transferred to secondary 
loop, or 



M X C C -^— = W c C c [T jx (t)-T ox (t)] 



(20) 



-U x A x [T x (t)-T 8 ] 



This equation was linearized by the same method 
used to obtain Equation 10; therefore, the connect- 
ing relationship 



t 



n o (t) = T c K c * (t) + K c / €(t) dt (23) 



and 



e(t)=T ref- T ove (t 



(24) 



The control system is made up of a comparator, a 
controller, and a control rod drive unit. The com- 
parator, defined mathematically by Equation 24, 
compares a reference temperature to an average 
temperature to form an error signal. The reference 
temperature is a predetermined function of demand- 
power level. Average system temperature in this 
study 



T ox<t) + T oc (t) 



(25) 



is an arithmetic average of high and low system 
temperature. If Equations 11, 21, 24, and 25 are 
combined, the error signal equation becomes 
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T ix (t) T jc (t) 

H + 



(26) 
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(21) 



is valid in describing the operation of the steam 
generator. A more convenient form of Equation 20 
is 



M*C ( 



dT x (t ) 
dt 



= 2W C C c T ix (t)-U X A X T S 
-K A x + 2W c C c ]T x (t) 



(22) 



the most convenient form of the error equation to 
use in the simulation. The controller is defined 
by Equation 23 and is self-explanatory. 

The control rod drive unit, which positions the 
control rods, can be defined by the transfer function 



K' 



S( l+r m S) 



(27) 



7oThe Control System — Automatic control of 
power level is achieved by positioning the control 
rods in accordance with the error signal repre- 
senting the difference between desired and actual 
power levels. Actual power output is proportional 
to neutron density and the power-demand signal 



The unusual forcing function, (n 0 -n)/n, is used to 
compensate in part for the non-linearities of the 
reactor. In practice, the reactivity contribution of 
control rod position is a non-linear function. In the 
interest of simplicity, it was assumed proportional 
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to control rod position. Therefore, Equation 27 
becomes 



n Q -n 



c s(|+T m s) 



(28) 



where K m is a modified gain constant. The most 
convenient representation of Equation 28 for simu- 
lation purposes is in terms of variable as de- 
fined by Equation 7. By combining Equations 7 and 
28 and converting the transfer function to a dif- 
ferential equation, the governing differential equa- 
tion for control rod drive unit becomes 



d 2 M t) 
dt 2 



I dMt) 
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dt 
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" n 0 (t)- n(t) 1 
n(t) J 



(29) 



The mathematical model for the entire system, 
which is summarized in Figure 1, has now been 
developed with the exception of the scram system. 
This system, which must be considered as manually 
controlled because of equipment limitations, is 
discussed in the mechanization section. No mathe- 
matical equations can be derived to define this 
system. Note that reactivity is highly negative 
when the scram rods are inserted in the reactor. 
Since these rods are spring-loaded or magnetically 
held, no time lag is experienced when they are 
inserted. 



MECHANIZATION OF MATHEMATICAL MODEL t 

The mathematical model was magnitude- scaled for 
a ±10 volt, PACE® TR-10 General Purpose Analog 
Computer, in the conventional manner using the 
data contained in Tables I and II. Table I is a sum- 
mary of the parameters encountered in this type 
of system. Table II is a summary of computer 
variables, initial values, and scale factors. 

TABLE I. SUMMARY OF PARAMETERS 



Fuel Elements, Etc. 
n m AHf / UA= I000°F 
n m AH f / M f Cf = 200°F/sec 
M f C f /UA = 5sec 

Plenum Chambers 



Mj /W c =2sec 



Steam Generator 
U X A X /U x A xt 2W C C C = 0.4 

U x A x /W c Cc= 4/3 

M x C c /U x A x +2W C C c = lsec 



350°F 



Reactivity 
cc= -2 x I0" 5 °F H 
Sk c = -0.005 (MAXIMUM) 



Reactor Core 

UA/W C C c =0.4 

M C /UA+2W C = l/3sec 
c c c c 

UA/UA +2W C C C = I /6 



Transport Delays 
^o/vo = 8 sec 
l\/ v j = 8 sec 



Control System 
T ref = I000°F (MAXIMUM) 
K* = I0 -4 (°F-sec) _l 
r r - 10 sec 

K m =0.05 sfic" 1 



Reactor Kinetics 
(U 235 Fuel) 

1 = 10 sec 



TABLE II. SUMMARY OF VARIABLES 



The quantities n 0 * (t), n m , K*, and n* (t), which 


Variable 


Initial Value 


Maximum Value 


Scale Factor 


Scaled Voltag 


appear in the remainder of this paper, are the 


T f 


1000° F 


2000°F 


1/200 


[T f / 200 ] 


result of redefining the neutron density as 


T c 


500°F 


I000°F 


1 / 100 


[t c / 100 ] 




T oc 


600°F 


1000° F 


1/100 


[t oc /ioo ] 


n( t) 


T,x 


600°F 


I000°F 


1/ 100 


[T jx / 100 ] 


n*(t)= (30) 


Tx 


500°F 


1000° F 


1/ 100 


[t x / 100 ] 


n m 


Tox 


400°F 


I000°F 


1/ 100 


[Tox /'GO ] 




To 


400°F- 


1000 °F 


1 / 100 


[t 0 / 100 ] 


where n* is reduced neutron density and n m is the 
maximum practical neutron density. Neutron den- 


Tic 


400°F 


I000°F 


1 / 100 


[T| c / IOO] 


n* 


0.5 


1.0 


10 


[lOn* ] 


sity was normalized to prevent scaling difficulties 


n*o 


0.5 


1.0 


10 


[lOn* 0 ] 


which could be encountered in mechanizing the 


n$-n* 


0 


0.5 


20 


[20(n5-n*" 


mathematical model. This also aids in presenting 


€ 


0 


200°F 


1/100 


[ * /IOO ] 


the results of the simulation in more general 


8k 


0 


0.005 


2000 


[2000 8k] 


terms. 










[2000 


H- 


0 


0.005 


2000 



tTTie subject of analog computer programming and scaling will not be 
discussed in this paper since several excel lent references are avail- 
able (3, 5). 
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Scaled voltage equations, with the exception of re- 
actor kinetics and transport delay equations, are 
summarized in Table III. 

The equations require no explanation. It should be 
noted that none of these equations required time 
scaling. The system was simulated in 4 'real time", 
therefore. 

TABLE III. SCALED COMPUTER EQUATIONS 



were also mounted. This RKN simulated the six 
delayed-neutron-equations, using six ' 'RC" passive 
element networks in the feedback loop of a high 
gain amplifier. When the RKN and the high gain 
amplifier are used in conjunction with a feedback 
capacitor, Cq, and an input resistor, Ro, Equa- 
tions 1 and 2 are completely simulated. This cir- 
cuit is shown in Figure 2. It accepts an input volt- 
age, which is the negative product of the voltages 



Control System 
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Reactivity 

[2000?k] = IO(20k)[lo]-l0^4 xlO 4 M^^oj + [2000 u] 
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Fuel Element, 
Coolant Heat Transfer 
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dt [\00\ V M cCj L 200 J \5M C C C / L'OOJ \ Mc C c j [ lOOj 
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Heat Transfer 
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dt boo J " \m 0 c c J L100J \M 0 c c y L100 J 

d_ \Jjc_] / Wc Cc \ [ Tj I / W c C c \ r T ic 1 
dt LlOO J \Mj C c / L 100 J I Mj C C J L 100 J 

Steam Generator 

d [ T x 1 / 2W C C c \ [ T ix 1 / U X A x T S V / U x A x -h 2W c C c \ f T x 1 
dt L 100 J " Ux C c / L 100 J \I000 M x cj \ M X C C / L 100 J 



- e 8k < 
10 



— VW- 



RKN 



Figure 2. Reactor Kinetics Circuit 



representing neutron density and reactivity, and 
has as its output the voltage representing the neu- 
tron density. It should be noted that seven ampli- 
fiers were saved by using the RKN. The values of 
the mean neutron lifetime and the maximum value 
of the reactivity determine the magnitude of Rq 
and Co. The theory of the RKN is discussed in 
Appendix B. 



NOTE: Quantities in [ ] = Scaled Variable, ( ) = Pot Setting 

1/2, 1, 10, etc. = Input gain when appearing on left of pot setting, reference 
voltage when appearing on right of pot setting 

The conventional PACE TR-10 General Purpose 
Analog Computer includes twenty amplifiers, 
twenty-four attenuators, and a variety of non- 
linear equipment. In order to simulate the reactor 
kinetics and transport delays some of this equip- 
ment was exchanged for special purpose compon- 
ents specially developed for nuclear studies. The 
four wirewound, 10-turn attenuators and two func- 
tion switches, normally located on the control 
panel, were removed, and replaced by a control 
panel for a Reactor Kinetics Network (RKN). This 
network was located in the rear of the computer. 
The control panel contained three two-position func- 
tion switches (one labeled " scram"), and their 
patching terminations. An initial- condition potenti- 
ometer and input and output patching terminations 



Representation of t ran sport del ay was accomplished 
by storing discrete samples of a time- varying volt- 
age in a capacitor memory, and subsequent read- 
out of the stored voltages to generate a delayed 
"staircase" approximation of the input function. 
The switching and storing functions of the memory 
were performed by four components, which must 
be added to the conventional TR-10 Computer. 
These components are: a timing unit, a ring counter, 
capacitor memories, and multiple sample-hold 
units. The timing unit, which replaces the repetitive 
operation drive unit on the TR-10 control panel, 
generates a pulse train to the ring counter. The 
ring counter, located in the nonlinear row of the 
computer, drives the capacitor memory and multi- 
ple sample hold units. These are also located in 
the non- linear row. Each transport delay requires 
a high gain amplifier, a capacitor memory, and a 
multiple sample hold unit. The transport delay 
circuit is shown in Figure 3. 
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j = n 



e Q (t) 



ej (t-r) 
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various reference temperature inputs. These runs 
were performed to make a steady- state check of 
the system. The system was designed for a 500°F 
reference temperature input to be compatible with 
the initial conditions and potentiometer settings. 
Therefore, if a reference temperature of 500°F 
was patched into the circuit, it should draw a hori- 
zontal line along the initial value of n*(5), which 
was one-half, or five volts. The steady- state check 
was successful as shown in Figure 5. The slight 
off- set in the steady state value is due to the ac- 
cumulated error of the computer components, each 
of which is individually rated at 0.10%. This off-set 
is about 1.0%. t 



Figure 3. Transport Delay Circuit 

The reactor kinetics and transport delay circuits 
designed for the TR-10 computer also are available 
on the TR-20 unit so that the program can be run 
equally well on either machine. 



Figure 5 shows that maximum reference tempera- 
ture (or power level) of the system lies between 
600°F and 700°F. During these runs, function 
switches in the computer circuit were in the "up" 
position, which meant that the control system was 
being used. 



The capacitor memory contains 16 storage capaci- 
tors, and relays associated with each capacitor. 
This unit, with the aid of the multiple sample-hold 
circuit and the timing equipment, is capable of 
sampling rates of 0.1 to 10 samples per second. 
Sampling rate is determined by a multiposition 
switch and vernier, located on the front of the 
timing unit. In the simulation under consideration, 
the time delav was 8 seconds: therefore n Rarnnlinp; 
rate of 2 samples per second was required. Trans- 
port delay units are wired to begin sampling simul- 
taneously with computer start-up. Provision is also 
made in this unit for obtaining an initial condition 
voltage at its output; however, the same voltage is 
stored on all sixteen capacitors. 

Please note that when a fourth order Pade approxi- 
mation is used to simulate a transport delay, it 
requires six amplifiers and six potentiometers. 

The computer diagram for the simulation of the 
primary loop is shown in Figure 4. The potentiom- 
eter settings and amplifier outputs associated with 
this circuit are tabulated in Table IV. This table 
consists of the static check potentiometer settings 
and amplifier outputs, as well as the 1 'steady- state 
check" potentiometer settings. The purpose of 
the steady- state check will be discussed in the 
results. 

RESULTS 



The same circuit was used also to obtain repre- 
sentative variations of n*(t),5 K(t), e(t), T f (t), T c (t) 
and T x (t) versus time (for a reference temperature 
of 650°F). Typical results are shown in Figures 6, 
7, and 8. A hazard evaluation of the control system 
is shown in Figure 6. This was accomplished by 
throwing the function switch feeding the final inte- 
grator of the control circuit (amplifier #18) to the 
down position after ii seconds of operation. This 
grounded the integrator- input and fixed the position 
of the control rods, which is the output of amplifier 
#18. The system stabilized itself by virtue of the 
temperature feedback of the fuel, but not at the 
desired level. The noise on the e(t) versus time 
curve, shown on Figure 7, is due to the switching 
in of the transport delay circuits as well as the 
quantitizing of the input function. This presents no 
problem if the output of the transport delay is fed 
to an integrator. The integrator, which acts as a 
smoothing device, produces a smooth curve. This 
can be proven by referring to the graph of T x (t) 
versus time. 

The third set of results, illustrated in Figures 9, 
10, and 11, show the behavior of the system for a 
step change in reactivity( T ref = 50 0°F), for manual 
control of the reactor, and for the effect of scram- 
ming, respectively. In the first case, coefficient 
potentiometer #1, which feeds amplifier #4, was 
set so that the voltage at the output of the amplifier 



The first series of computer runs, performed after 
a successful static check, resulted in a family of 
curves of reduced neutron density versus time, for 



t Although the authors have not performed this simulation on the 
TR-20 computer as of this writing, it is likely that the greater com- 
ponent accuracy of the TR-20 machine would resu It in a substantial ly 
lower total error. 
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was -5 volts, which corresponds to a step change 
in reactivity of -2.5 x 10~ 3 . As shown in Figure 9, 
the control system returned the neutron density 
to the initial value rapidly and with typical response. 

Manual control of the reactor was illustrated by 
removing the control system from the computer 
circuit. This was accomplished by putting both 
function switches in the "down" position. The com- 
puter then produced the initial value of n* (t) illus- 
trated by the ' { steady- state check", and by the 
first 30 seconds of the plot shown in Figure 10. 
At this point, coefficient potentiometer #1 was 
varied manually so as to increase the power level 
of the reactor from a neutron density level of 0.50 
to 0.80. 

A situation calling for the use of the scram rods, 
which are assumed to be present in addition to the 
control rods, was simulated by setting the reference 
temperature to an arbitrary value greater than 
700°F, As previously shown, 700°F is too large for 
normal reactor behavior; therefore, the flux level 
proceeded to a maximum level as shown in Figure 
11. When the flux reached its maximum value, the 
two function switches and the " scram'' switch 
(which feeds the reactivity summer — amplifier 



#4), were thrown simultaneously. This removed 
the control system from the circuit, and forced 
the reactivity to its lower limit of -10 volts. This 
was accomplished using a diode limiter on ampli- 
fier #4, as shown in Figure 4, and placing a large 
voltage on its input via the scram switch. The 
neutron level responded as expected by dropping 
quickly toward zero, first rapidly and then grad- 
ually. The same procedure was repeated for the 

T (t)andT (t) versus time, and their behavior was 

x w c 
as expected. 

CONCLUSIONS 

The results, as shown, are sufficiently accurate 
for an approximate analysis of over- all system 
behavior and are of use in practice. 

This simulation was possible on a twenty- amplifier 
analog computer only through the use of the special 
Reactor Kinetics Network and Transport Delay 
Simulator components. Their use reduced the re- 
quired number of amplifiers and associated po- 
tentiometers by at least seventeen, and provided 
the necessary special circuitry for reactor kinetics 
and time delay simulations. 
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TABLE IV 

TR-10 POTENTIOMETER ASSIGNMENT SHEET 



PROBLEM NUCLEAR REACTOR 



PHT 

rU 1 

NO. 


DADAMCTCD 

DESCRIPTION 


SETTING 
STATIC 
CHECK 


SETTING 
RUN NO.l 
* 


SETTING 

RUN 
NUMBER 1 


SETTING 
RUN NO. 


NOTES 


POT 
NO. 


1 


20 k 


0.400 


0.400 






TRIM POT 


1 


2 


T REF /\000 


0.550 


0,500 








2 


3 


2W C C C /M X C C 


0.600 


0.600 








3 


4 


LIMITER 


0.500 


0.500 






LIMIT TO-IOV 


4 


5 


UxAxT s /l0 3 M x C c 


0.140 


0.140 








5 


6 


T x C0)/|000 


O.800 


0.500 








6 


7 


T K (0)/\000 


0.600 


OAOO 








7 


8 


2 n * (0)/lO 


O.300 


0.5 OO 








8 


9 


SCALING RATIO 


O.2O0 


O.ZOO 








9 


10 


W C C C /5M C C C 


0.2 50 


0.Z50 








10 


II 


UA/M F C F 


O.ZOO 


O.ZOO 








1 1 


12 


(UA + 2W c Cc)/IOM c C c 


0.50O 


O.30O 








12 


13 


T F (0)/2000 


0.6 00 


O.50O 








13 


14 


4 X.I0 4 /^| 


0.800 


0.800 








14 


15 


T 0 (0)/!000 


0.600 


0.600 








15 


16 


K m /Tm 


O./OO 


O.lOO 








16 


17 


i/iox m 


0.200 


O.200 








17 


18 


T C (0)/IOOO 


0.800 


0.500 








18 


19 


ZOO K"T C 


O.ZOO 


O.ZOO 








19 


20 


2000 K * 


O.ZOO 


O.ZOO 








20 


21 


n*(o) 


0.500 


0.500 






RKN IC POT 


21 


22 














22 


23 














23 


24 














24 
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TABLE V 

TR-10 AMPLIFIER ASSIGNMENT SHEET 



PPORI FM NUCLEAR REACTOR 



AMD 

AJVIr 
NO. 


FB 


ni ITDI IT 
UU 1 rU 1 

VARIABLE 


STATIC CHECK 


NOTES 


CALCULATED 


MEASURED 


CHECK PT. 


OUTPUT 


CHECK PT. 


OUTPUT 


1 


I 


-io" 3 fedt 


-0.45 


-10.00 


✓ 


✓ 




2 


r 


- [Tx//00] 


-3.0 


-8.0 


✓ 


y/ 




3 


s 


- [lOn*] 




-5.0 




✓ 




4 


s 


+ [zooo b k] 




-8.0 




✓ 




5 


5 


- [e/ioo] 




+ 4.5 




y 




6 


5 


- [ZOOO S kj 




+ 8.0 




v/ 




7 


S 


+ [e/ioo] 




-4.5 




y 




8 


5 


+ [2(n<*-n*)] 




-0.9 








9 


HG 


+ [Txx /IOO] 




+ 6.0 




✓ 




10 


HG 


-i [Tx/100] 




-10/3 




y 




II 


I 


+ [T F /I00] 


-o.i 


+ 6.0 


y 


y 




12 


I 


-[Tc/IOO] 


-3.0 


-8.0 


/ 


y 




13 


I 


+ [T o /l00] 


-2.0 


+ G.O 


✓ 


y 




14 


I 


+ [Tic/lOOj 


-2.0 


+ 6.0 


✓ 


y 




15 


HG 


L TJX" J 




+ 1.8 




y 




16 


5 


*[^\:""] 




- 1.8 




y 




17 


I 




+21.8 


+ io 


✓ 


y 


-|OV TO I.C. IN 
STATIC CHECK 


18 


I 












USED AS CHECK 
AMPL. DURING S.C. 


19 


HG 


-20O0 ^*Sk 




+ 4.0 




y 




20 


HG 


+ [10**] 




+ 5.0 




y 
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Figure 5. n*(t) Versus t for Various Reference Temperature (T Q ) Inputs 
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Figure 6. The Effect of Control System Failure on Reactor Operation 
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t in seconds 



Figure 7. e(t) and 8k(\) Versus t 



in °l 




t in seconds 



Figure 8. T c (t), T x (t), and T f (t) Versus t 




f in seconds 

Figure 9. The Effect of a Negative Step Input in Reactivity on Reactor Operation 
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t in seconds 

Figure 10. n* (t) Versus t for Manual Control of the Reactor 
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Figure 11. The Effect of "Scramming" on n* (t), Tf (t), and T x (t) 
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APPENDIX A 



LINEARIZATION OF FLUID HEAT TRANSFER EQUATION 



Consider a fluid flowing through a closed conduit which is gaining thermal energy from its surroundings. If it is assumed 
that the ambient temperature of the surroundings is not a function of conduit length, the temperature history of the fluid 
may be represented by the familiar equation 

AVc ^ar^t) =u v a [ Ta(t) _ T(M) ]_ wc (A-i) 

This equation is governed by the initial conditions 

T(o,x) = f (x) (A-2) 

and 

T(t,o) = g(t) (A-3) 
If one assumes that a linear temperature profile with respect to distance, 

T(x,t) = T(t,o) + (T(t,L) -T(t,o)) 1 (A-4) 
exists at any time, t, the mean fluid temperature^f (t), is 

- (f)= ^T(x.t) dx T(t,o) +T(t,L) (A _ 5) 



dx 



2 



This equation relates the inlet, outlet, and mean temperatures of the fluid which are all time varying functions. 
Integrating Equation A-l with respect to distance yields 

ApC jl axix^ dx = Ua ; o L [ Ta(t) _ T(xt) ] dx _ wc / o L iI£i! dx (a-6) 

or 

LApC / 0 ' iV" [ T ^> o) + (T(t,L)»T(t l o) > ) s] ds = UaL [t q (T) [T(t,o)+ (l(t,L) - T (t,o) )] xds 

- WC [t (t,L) -T(t,o) ] 

which reduces to 

LA/oC $5±& = UaL [t q ( t) — T( t)] — WC [T(t,D — T(t,o)] (A -8) 

Note that 

LA p = M (A-9) 

and (A-10) 

□ L = A h 

where M is the mass of the fluid in the conduit and A^ is the heat transfer area. Substituting these expressions in Equa- 
tion A-8, we obtain 

MC dllL) = UA h [ T Q (t)-T (t) ] -WC [ T(t,L) - T(t,o)] (A-ll) 
which, except for notation, is Equation 10. 
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APPENDIX B 



REACTOR KINETICS NETWORK THEORY 

The passive network for simulating the reactor kinetic equations is shown in simplified form in Figure B-l. If the input 
stage of the amplifier draws negligible current*Kirchoff's current law gives 

6 

i + i f + I I. =0 (B-l) 

If one further assumes that the summing junction is at zero potential, Equation B-l can be written 

6 

o-TT- i- 2ii= 0 (B-2) 



Ro +Co H-^, li=0 



If the voltage across capacitor C. is called e., and its charge Q., then with sign conventions as shown in Figure B-l, we 
obtain 3 3 3 

Qj =Cj ej 

(B-3) 

ej + ij Rj = e 0 (B-4) 
dQj (B-5) 

'i = "dT 

In deriving B-4, we again used the assumption that the summing junction is at zero potential. 



Eliminating e. from B-3 and B-4: 



e, 



b _ W J 
J " Rj Rj Cj 



Q : (B-6) 



Substituting B-6 into B-2: 



where 



Solving B- 7 for 



ej de 0 | 6 Qj 

* 0 +c °ir + ^ e °-}^r° (B " 7) 

R * Rj 

de 0 
dt 

de o ~ e i _ J_ J_ | _Qj_ 

dt " R 0 C 0 RC 0 60 C 0 f Rj Cj < B -8) 



Eliminating i from B-5 and B-6, we obtain 



d Qj Qj 

dT + R-C- = Rj" 6 ° (B " 9 > 



Comparing B-8 and B-9 with the reactor kinetic equations 



dn n 8k p v X- r • (B-10) 

JT = IT " T* n + S X J C i (B } 

dC J x „ P \ (B-ll) 



dt 
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it becomes clear that if we make (-ei) proportional to n5 k, and e Q proportional to n, the equations will have the same 
form. The values of Rj, Cj, R 0 , and C Q can be calculated in terms of the known constants Xj, 0j, and i*, and the 
scale factors chosen for n and«5k. 

It turns out that the Rj and Cj values may be fixed if a given fuel is used, and C 0 and R 0 chosen as functions of / * (which 
may vary widely with reactor design) and the scale factor for 5k (which depends on the estimated maximum value of 5k). 
The six resistors and capacitors, Rj and Cj.are packaged within the reactor kinetics unit, while R Q and C 0 are patched 
externally, so that various values may be used to suit a particular problem. 

Typical values of I * and <5k m (the maximum reactivity for the simulation) are given in Table 1-B, along with the calcu- 
lated values of R 0 and C Q . The component values shown in this table are available from Electronic Associates, Inc. as 
plug-in units. 

R 6 + 
— WV 11 



TABLE 1-B. SELECTED VALUES OF R AND C t 

o o 

ak^xio 2 



I* 

seconds 

io- 3 

ID" 4 

io" 5 



10 



10 



10 



10 



-4 



10 



-4 



10 



m 

dimensionless 
0.05 
0.05 
0.05 
0.20 
0.20 
0.20 
0.50 
0.50 
0.50 



o 

megohms 
10.0 
10.0 
10.0 
2.5 
2.5 
2.5 
1.0 
1.0 
1.0 



^o 



mfd 

0.200 

0.020 

0.002 

0.200 

0.020 

0.002 

0.200 

0.020 

0.002 



R o 



"6 



R 5 

-WV— 



" 4 

-AAA/— 



c 5 

, + 



14 



R 3 
-A/vV- 



'2 



R , 

A/W- 



SJ 



t Electronic Associates, Inc.; PIR 6239 



Figure B-l. Simplified Diagram of the 
Reactor Kinetics Network 
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APPENDIX C 



SYMBOLS FOR COMPUTING CIRCUITS t 



POTENTIOMETERS 



GROUNDED 



UNGROUNDED 



-0- 



o < k < 1 



e 0 = e, + k(e 2 -e, ) 



OPERATIONAL AMPLIFIERS 



INTEGRATOR 



10 




e o =e ic-l < e i +10e 2 )dt 



SUMMER 



0.1 



6 0 = ~^ e ' + °' le 2 + ,0e 3^ 



HIGH -GAIN 



/i > > o 



e 0 =-/a ej 



MULTIPLIERS 



x 1 s 



(vn 



-XY 
°" 10 



FUNCTION SWITCH 



-e, e 0 =e, U- POSITION 

62 e o =e 2 D-POSITION 



t TR-10 Operators Manual 
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APPENDIX D 



TABULATION OF COMPUTER COMPONENTS 

The following TR-10 Computing Components (which were housed in a standard TR-10 Console), patching 
sories, and readout devices were required to carry out the simulation. 

Computing Components 



20 - 


Operational Amplifiers 


8 - 


Integrator Networks 


2 - 


Quarter- Square Multipliers 


2 - 


Capacitor Memory Units 


2 - 


Multiple Sample Hold Units 


1 - 


Ring Counter 


1 - 


Timing Unit 


20 - 


Coefficient Setting Potentiometers 


1 - 


Reactor Kinetics Network and Control Panel 



Patching Components 
122 - Patch Cords 

22 - Black, 6 inch 
67 - Brown, 12 inch 
30 - Orange* 18 inch 

3 - Blue, 30 inch 

4 - Multiple Connectors 
11 - Feedback Resistors 

2 - 100K 

9 - 10K 

38 - Input Resistors 

4-1 megohm, includes input resistor to RKN 

4 - 200K 
7 - 100K 

27 - 10K 

1 - Diode 

1 —Capacitor, 0. 02 fd feedback for RKN 
Readout Devices 

1 - X-Y Plotter (with time base generator) 
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LIST OF SYMBOLS 



MEANING 
Heat transfer area in the reactor 
Heat transfer area in the steam generator 
Cross sectional flow area 
Specific heat of the coolant 
Specific heat of the fuel, moderator, etc. 
Concentration of neutrons in delayed neutron group "j" 
Fluid specific heat 
Reactivity constant 
Control rod drive unit gain 

Controller gain 

Mass of coolant in the reactor 
Mass of fuel, moderator, etc. 
Mass of coolant in the inlet plenum chamber 
Mass of coolant in the outlet plenum chamber 
Mass of coolant in the steam generator 
Operator 
Ambient temperature 

Average coolant temperature in the reactor 
Average fuel temperature 

Coolant temperature at the input to the reactor inlet 
plenum chamber 

Coolant temperature at the outlet of the reactor outlet 
plenum chamber, or temperature coefficient of reactivity 
reference temperature (temperature at which temperature 
contribution is zero). 

Average temperature of secondary fluid in steam 
generator 

Average coolant temperature in the steam generator 
Temperature of coolant entering the reactor 
Temperature of coolant entering the steam generator 
Temperature of coolant leaving the reactor 
Temperature of coolant leaving the steam generator 
Reference temperature 
Average system temperature 
Fluid temperature 

Overall coefficient of heat transfer in the reactor core, 



UNITS 
Sq Ft 
Sq Ft 



Neu/sec 

Dimensionless 

Ft/sec 

Neutrons 
•F- sec-cm 2 

Lbs 

Lbs 

Lbs 

Lbs 

Lbs 

Seconds * 



SYMBOL 
U' 
V 

w 

c 

w 

X 



li 

lo 
I* 



5k 

c 

5 k, (o) 

5k f 

5k 



e(t) 
M(t) 



Control rod reactivity variable 



Overall coefficient of heat transfer in the steam 
generator 



sec -Fr 
Dimensionless 

Dimensionless 



MEANING UNITS 

Overall coefficient of heat transfer Dimensionless 

Mean fluid velocity Ft/sec 

Mass flow rate of coolant Lbs/sec 

VA, mass rate of flow of the fluid ■■ 

Control rod position Dimensionless 

Heat transfer area per unit length of conduit 

Index for delayed neutron groups Dimensionless 

Effective length of inlet piping system Ft 

Effective length of outlet piping system Ft 

Effective neutron lifetime Seconds 

Neutron density Neutron/cucm 

Reduced neutron density Dimensionless 

Maximum practical neutron density Neutron/cucm 

Demand-power level neutron density Neutron/cucm 

Time Seconds 

Mean velocity of coolant in outlet piping system Ft/ sec 

Mean velocity of coolant in inlet piping system Ft/ sec 

Position along conduit 

Temperature coefficient of reactivity °F-1 
Fraction of prompt neutrons appearing in delayed 

neutron group "j" Dimensionless 

Reactivity Dimensionless 

Reactivity contribution of control rod positions Dimensionless 

Initial reactivity contribution of control rods Dimensionless 

Built-in reactivity of fuel Dimensionless 

Reactivity contribution of reactor poisons Dimensionless 

Reactivity contribution due to the fuel temperature Dimensionless 

Heat of fission Btu - cu cm 

Sec - Neu 

Error signal °F 

Decay constant associated with group "j" Dimensionless 
Departure of control rod reactivity from its initial value Dimensionless 

Outlet piping system delay time Seconds 

Controller time constant Seconds 

Inlet piping system delay time Seconds 

Control rod drive unit time constant Seconds 

Fluid density 
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